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F I G. 2 



t i 


=Tc i +Tm i 




(1a) 


t i 


: D 1 FFUS 1 ON-CODE TRANSM 1 SS 1 ON 


TIME 




Tc 


i : VALUE OF D 1 G 1 TS EXPRESS 1 NG 


A NUMBER 


SMALLER 




THAN 1 msec 








(PHASE OF DIFFUSION CODE) 






Tm 


i : VALUE OF D I G I TS EXPRESS I NG 


A NUMBER 


EQUAL 




TO OR GREATER THAN 1 msec 






ti 


=Tci+Ti+to 




(1b) 


ti 


: D I FFUS I ON-CODE TRANSM I SS I ON 


TIME 




Tc 


i : VALUE OF D I G I TS EXPRESS I NG 


A NUMBER 


SMALLER 




THAN 1 msec (PHASE OF DIFFUSION CODE) 




(T 


i +to) : VALUE OF D I G I TS EXPRESS I NG A NUMBER EQUAL 




TO OR GREATER THAN 1 msec 




to 


: REFERENCE TIME (COMMON TO SATELLITES) 




Ti 


: D I FFERENT I AL TIME 
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F I G. 3 



T 1 VAU, 




. to) = 


I Xi (ti)-Xo | -c { (tr- 
| X i (Tc i +T i +to) -Xo | 
-c {tr-(Tci+Ti+to)-»- 


r } 


s 

=0 


X i (t i ) 


= (x i (t i ) y i (t i ) z i (t i ) ) T 




■-■(3) 


AO — ^XO 


yo zo)^ 








f 1 (Xo. 


r 


. to) = 


| X1 (Tc1+T1+to)-Xo | 
-c {tr-(Tc1+T1+to) + 


T } 


=0 j 
■■■(2a) 


f 2 (Xo. 


r 


. to) = 


| X2 (Tc2+T2+to) -Xo | 
-c {tr-(Tc2+T2+to)+ 


T } 


=0 

■••(2b) 


f 3 (Xo. 


r 


. to) = 


| X3 (Tc3+T3+to) -Xo | 
-c {tr- (Tc3+T3+to) + 


r } 


=0 

■--(2c) 


f4 (Xo. 


r 


, to) = 


| X4 (Tc4+T4+to) -Xo | 
-c {tr-(Tc4+T4+to)+ 


r } 


=0 

•--(2d) 


f 5 (Xo, 


r 


. to) = 


| X5 (Tc5+T5+to) -Xo | 
-c {tr-(Tc5+T5+to)+ 


T } 


=0 

••-(2e) 
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F I G. 4 



POINT OF INTEREST 



(Xo(n) , r (n),to(n)) (11) 



A f 1 = 3 9 xo'(n) A xo+ 3yo'(n) A yo+ 9 a Z o'<n) Azo 



f i (Xo (n) . r (n) , to (n) ) + A f i =0 (13) 

(AXo, A t , A to) (14) 

AXo=(Axo Ayo Azo)T (15) 



9 f i - 9 f i 
9 xo (n) 9 xo 



xo=xo (n) (16) 



F(Xo(n), r (n) , to(n))+AF 
=F (Xo (n) , r (n) , to (n) ) +J (Xo (n) . v (n) . to (n) ) A P 
=0 (17) 



F=(f1 f2 f3 f4 f5) T (18) 

AP=(Axo Ayo A r Ato) T (19) 
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F I G. 5 



J (Xo (n) , z (n) , to (n) ) = 



afi 
3 xo (n) 

a xo (n) 

af3 

a xo (n) 

af4 

a xo(n) 

af5 



afi 

a yo (n) 

3f2 
a yo (n) 

3 yo (n) 

3f4 

a yo (n) 

af5 



afi afi afi 

a zo (n) a r (n) a to(n) 

af2 a f2 a f2 
a zo(n) a r (n) a to(n) 

a f3 af3 a f3 

a zo(n) a r (n) a to(n) 

a f4 af4 a f4 

a zo(n) a r (n) a to(n) 
a f 5 a f 5 a f 5 



afi 


xi (Tci+Ti+to(n)) 


-xo (n) 


a xo (n) 


I Xi (Tci+Ti+to(n)) 


-Xo(n) | 


afi 


yi (Tc i +T i +to (n) ) 


-yo (n) 


ayo(n) 


I Xi (Tci+Ti+to(n)) 


-Xo(n) | 


afi 


zi (Tc i +T i +to (n) ) 


-zo (n) 


a zo (n) 


I Xi (Tci+Ti+to(n)) 


-Xo(n) | 



afi 

a ro(n)" 0 



$ xo (n) a yo (n) a zo (n) a z (n) a to (n)„ 

(21) 



(22) 
(23) 
(24) 
(25) 



afi _ (Xi (Tci+Ti+to(n))-Xo(n))-Vi (Tc i +T i +to (n) ) 
a To(n) - | Xi (Tci+Ti+to(n))-Xo(n) | +L 

A-B=(Xi (Tci+Ti+to(n))-Xo(n)) -Vi (Tci+Ti+to(n)) 



(26) 
(27) 
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F I G. 7 

NEXT POINT OF INTEREST 

(Xo(n+1), z (n+1) . to(n+D) 

= (Xo(n)+AXo, z (n) + A r,to(n)+Ato) 

s (33) 

CONVERGENCE CR I TER I A ( £ : CONSTANT) 

I AXo | <£ (35) 

I AP |<e (36) 

| F (Xo (n) . r (n) . to (n) | < e (37) 
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F I G. 8 

CALCULATION OF DIFFERENTIAL TIME Ti 

t1=Tc1+to (1c) 

ti=Tci+Ti+to (i^2) (1d) 

T i = (T i +to) -to= (t i -Tc I ) - (t1 -Tc1 ) (41 ) 

c(tr-ti)= | Xi (ti)-Xo | (42) 

tj=tr _ |Xi(ti)-Xo| (43) 

t1=tr _ j X1 (tl)-Xol (44) 

C 

T i = (i X1 (t1)-Xo| +Tc1 )_( |Xi(ti)-Xo| +Tc .) 

(45) 

cx1.0x10" 3 [m]~3.0x10 5 [m] (46) 

ti=t 1-2. 0 x 10~ 2 [SECONDS] (47) 
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F I G. 9 



TOLERANCE OF APPROXIMATE POSITION OF RECEIVER 



cTi = ( | X1 (t1)-Xo | +cTc1) 

-( I Xi (ti)-Xo | +cTci) (51) 



( | X1 (t1)-(Xo+AX) | +cTc1) 

-( I XI (ti)-(Xo+AX) | +cTci) 

= ( | X1 (t1)-Xo | + | AX I +cTc1) 

- ( | X i (t i ) -Xo | - | A X | +cTc i ) 

= ( | X1 (t1)-Xo | +cTc1) 

-( | Xi (ti)-Xo | +cTci)+2 | AX | 

=cTi+2 | AX | 



| AX | : DISTANCE BETWEEN COORDINATES 
REPRESENTING APPROXIMATE 
POSITION OF RECEIVER AND ACTUAL 



POSITION OF RECEIVER 
I A | > | B | (63) 

|A|-|B|^|A-B|^|A| + |B| 

(64) 

I Xi (ti)-(XcH-AX) | 

= I (Xi (ti)-Xo)-AX | (65) 

A=X i (t i ) -Xo (66) 

B=AX (67) 

I X i (t i ) -Xo | - | A X | 
^ I (X i (t i ) -Xo) - A X | 

^ | X i (t i ) -Xo | + |AX| (68) 
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F I G. 1 0 



TOLERANCE OF ERROR OF RECEIVER INTERNAL TIME 



t1 + At=ti + At+2.0x10" 2 [SECONDS] (53) 



( | X1 (t1 + At)-Xo | +cTd) 

-( | Xi (t1 + At)-Xo | +cTci) 
= ( | X1 (t1 + At)-Xo | +cTc1) 

- ( | X i (t i + A t+2. 0 x 1 0" 2 ) -Xo | +cTc i ) 
= ( | X1 (t1)-Xo | +At x 1 . 0 x 10 3 +cTc1) 

-( I Xi (ti)-Xo | - A t+2. Ox 10 -2 ) x 1. Ox 10 3 +cTci) 
= ( I X1 (t1)-Xo | +cTc1)-( | Xi (ti)-Xo | +cTci) 

+ (2 A t+2. 0x10 2 ) x1.0x10 3 
=cTi + (2At+2. Ox 10" 2 ) x 1. Ox 10 3 (54) 

A d= (2 A t+2. Ox 10 2 ) x 1. Ox 10 3 [m] 

~2Atx 1.0x10 3 [m] (55) 
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F I G. 1 1 

2 | AX | +2 | At | x 1. Ox 10 3 <1. 5x 10 5 (56) 
| AX|+1.0x10 3 x | At | <7. 5 x 10 4 (57) 
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F I G. 1 2 



SECOND METHOD 

f1 (Xo, r . to)= | X1 (Tc1+T1+to)-Xo | 

-c {tr-(Tc1+T1+to)+ r } =0 

(2a) 

f 2 (Xo, r , to) = | X2 (Tc2+T2+to) -Xo | 

-c {tr-(Tc2+T2+to)+ r } =0 

(2b) 

f3(Xo, r . to)= | X3 (Tc3+T3+to) -Xo | 

-c {tr-(Tc3+T3+to) + r } =0 

(2c) 

f 4 (Xo. t , to) = | X4 (Tc4+T4+to) -Xo | 

-c Itr-(Tc4+T4+to)+ r } =0 

(2d) 

f5(Xo, r,to)= | X5 (Tc5+T5+to) -Xo | 

-c {tr-(Tc5+T5+to)+ r } =0 

(2e) 

f6(Xo, r,to)= | X6 (Tc6+T6+to) -Xo | 

-c {tr-(Tc6+T6+to)+ r } =0 

- (2f ) 

AP=(J T J) 1 J T (-F) (61) 
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THIRD METHOD 

f 5 (Xo) = | Xo | - | Xoa | =0 (71 ) 

Xoa : APPROX I MATE POSITION OF RECEIVER 

d f 5 xo (n) 

a xo (n) | Xo (n) I 

3 f 5 _ yo (n) 
a yo (n) | Xo (n) | 

d f 5 zo (n) 
6 zo (n) | Xo (n) | 



(72) 
(73) 
(74) 
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F I G. 1 4 

C START ) 

_~3 — 



40 



41 



DETECT DIFFUSION- 
CODE PHASES OF 
SATELLITES WHOSE 
SIGNAL ARE BEING 
RECEIVED 






PERFORM 






POS I T I ON I NG 






CALCULAT I ON 






PROCESSING BASED 






ON CONVENTIONAL 






METHOD 




JL 44 



'SOLUT I ONJ 
JDBTA I NED? 



YES 



NO 



Sa:SET OF SATELLITES 

WHOSE SIGNALS ARE BEING 
RECI VED 

Sc:SET OF SATELLITES 
WHOSE DIFFUSION- 
CODE PHASES (Tc i ) HAVE 
BEEN CONFIRMED 

Sr:SET OF SATELLITES 
WHOSE DIFFUSION- 
CODE TRANSMISSION TIMES 
(ti) HAVE BEEN CONFIRMED 

SrcSccSa 




PERFORM 
POS I T I ON I NG 
CALCULAT I ON 
PROCESSING BASED 
ON THIRD METHOD 




PERFORM 
POSITIONING 
CALCULAT I ON 
PROCESSING BASED 
ON FIRST METHOD 




PERFORM 
POS I T I ON I NG 
CALCULAT I ON 
PROCESSING BASED 
ON SECOND METHOD 



YES 



49a 

f 



DETERMINE THAT 
POS I T I ON I NG 
WAS SUCCESSFUL 



( END ) 




DETERM I 
POS I T I ON 
WAS 



THAT 
ING 

UNSUCCESSFUL 
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50 



(start ) 



CALCULATE Ti 



k=1 

P=(Xoa. 0, toa) 



-80 



51 



FIND MATRIX 


J(P) BY 


USING Eqs. (22) TO (26) 




53 






FIND F(P) BY USING Eq. 


(2) 





52 



54 



FIND BY SOLVING 

Eq. (17) , THAT IS. 
-F (P) =J (P) Z^P 



YES 




DETERMINE THAT 
POS I T I ON I NG 

WAS SUCCESSFUL 



56 

^ 



k=k+1 




P=P+Z^P 



59b 



DETERMINE THAT 

POS I T I ON I NG 

WAS UNSUCCESSFUL 
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60 



(start ) 



CALCULATE Ti 



k=1 

P= (Xoa, 0, toa) 



80 



61 



62 



FIND MATRIX J (P) BY 
USING Eqs. (22) TO (26) 



63 



FIND F(P) BY USING Eq. 
(2) 



64 



FIND 2\P BY SOLVING 
Eq. (61) , THAT IS. 

^,P= (J (P) T J (P) )-1 J (P) T (-F (P) ) 




DETERMINE THAT 
POS I T I ON I NG 
WAS SUCCESSFUL 




DETERMINE THAT 

POS I T I ON I NG 

WAS UNSUCCESSFUL 



( END ) 
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( START ) 



70 



CALCULATE Ti 



k=1 

P=(Xoa, 0, toa) 



80 



71 



FIND MATRIX J (P) BY 
USING Eqs. (22) TO (26) 
AND (72) TO (76) 



72 



73 



FIND F(P) BY USING Eq. 
(2) AND (71) 



74 



FIND ^P BY SOLVING 
Eq. (17) , THAT IS, 
-F (P) =J (P) Z^P 



76 



k=k+1 




DETERMINE THAT 
POSITIONING 
WAS SUCCESSFUL 



DETERMINE THAT 

POS I T I ON I NG 

WAS UNSUCCESSFUL 
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( START ) 



F I ND X1 (toa) 



-81 



T1=0 

1X1 (toa)-Xoa | +Tc1 



-82 



i=2 



-83 



FIND Xi (toa) 



Tj=d1 _|X1 (toa)-Xoa| +Tc . 



80 




( END ) 
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F I G. 2 2 A 



n_F=\_Tl t i : D I FFUS I ON-CODE TRANSM I SS I ON T I ME 
ULJU X i (t i 



7 



i ) : POS I T I ON (3-D I MENS I ONAL COORD I NATES) 



Si 



C; VELOCITY OF LIGHT IN VACUUM 



tr : DIFFUSION-CODE RECEPTION TIME 
t : CLOCK ERROR 

Xo : POS I T I ON (3-D I MENS I ONAL COORD I NATES) 



F I G. 2 2 B 



I X i (t i ) -Xo | -c { (tr-t i) + r } =0 
Xi (ti) = (xi (ti) yi (tl) zi (ti)) T 
Xo=(xo yo zo)^ 



(91) 
(92) 
(93) 



I X1 (t1 ) -Xo | -c { (tr-t1 ) + t } =0 (91a) 

| X2 (t2) -Xo | -c { (tr-t2) + r } =0 (91b) 

I X3 (t3) -Xo | -c { (tr-t3) + r } =0 (91c) 

| X4 (t4) -Xo | -c { (tr-t4) + zr } =0 (91 d) 

f 4 (Xo) = | Xo | - 1 Xoa | =0 (95) 



Xoa : APPROX I MATE POSITION OF RECEIVER 



